The boson-fermion model: An exact diagonalization study 
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The main features of a generic boson-fermion scenario for electron pairing in a many-body corre- 
lated fermionic system are: i) a cross-over from a poor metal to an insulator and finally a supercon- 
ductor as the temperature decreases, ii) the build-up of a finite amplitude of local electron pairing 
below a certain temperature T* , followed by the onset of long-range phase correlations among elec- 
tron pairs below a second characteristic temperature T$, iii) the opening of a pseudogap in the DOS 
of the electrons below T* , rendering these electrons poorer and poorer quasi-particles as the tem- 
perature decreases, with the electron transport becoming ensured by electron pairs rather than by 
individual electrons. A number of these features have been so far obtained on the basis of different 
many-body techniques, all of which have their built-in shortcomings in the intermediate coupling 
regime, which is of interest here. In order to substantiate these features, we investigate them on 
the basis of an exact diagonalization study on rings up to eight sites. Particular emphasis has been 
put on the possibility of having persistent currents in mesoscopic rings tracking the change-over 
from single- to two-particle transport as the temperature decreases and the superconducting state 
is approached. 

PACS numbers: 74.20.-z, 74.78. Na, 74.20.Mn 



I. INTRODUCTION 

A great interest has recently been devoted in con- 
densed matter theory to models for interacting boson- 
fermion systems. Among the topics to which these mod- 
els apply, we recall the hole pairing in semiconductors 1 , 
the isospin singlet pairing in nuclear matted, d-wave hole 
and antifcrromagnetic triplet pairing in the positive-?/ 
Hubbard system^ (and possibly also in the t-J model), 
and entangled atoms in squeezed states in molecular 
Bose-Einstein condensation in traps 4 . 

In the present paper we consider a boson-fermion 
model (BFM) which captures a number of basic physi- 
cal properties of strongly interacting many-body systems 
where resonant pair states of bosonic nature form inside 
a reservoir of fermions. This BFM was originally con- 
ceived as a possible alternative and extension to the idea 
of bipolaronic superconductivity 5 - a scenario for systems 
with extremely strong electron-phonon coupling where all 
the electrons exist in form of locally bound pairs of small 
polarons. Since in such a case of extreme antiadiabaticity 
the mobility of bipolarons is expected to be vanishingly 
small, they have to be considered, for all intents and pur- 
poses, as remaining localized rather than condensing into 
a superfluid phase. An alternative to this situation was 
to consider a case where the electron-phonon coupling is 
of intermediate strength, such that the energy of locally 
bound isolated bipolarons is comparable with the gain in 
energy of itinerant electrons. This leads to a scenario 
where two-electron states fluctuate between pairs of itin- 
erant single-electron states close to the Fermi level and 
localized bipolaron states which pin this Fermi level. A 
detailed discussion on such a polaron-induced scenario 



for this BFM has been given recently 6 . Treating the 
bipolarons as bosonic entities which commute with the 
fermionic electrons, the Hamiltonian for such a system is 
given by 

H = (D - /,) + ( A b 2 M ) U + I) 

+ t c ^ c ^ +9 Y (W^JAT +Pl4\4i) •(!) 

The bipolarons are treated here as on-site entities hav- 
ing hard-core boson features and thus are represented 
by pseudo-spin operators {pf , p~ , pf}. cf a and de- 
note the creation and annihilation operators of the itin- 
erant electrons and g is the strength of the boson-fermion 
charge exchange interaction. The hopping integral for 
the itinerant electrons is given by t with a half band- 
width equal to D — zt, z denoting the coordination num- 
ber of the underlying lattice. The energy level of the 
hard-core bosons is denoted by Ab ■ A common chemical 
potential \i ensures the conservation of the total charge 
n-tot — fiF| + n Fl + ^riB {tib and n_F|,x denote the total 
number of the hard core-bosons and of the electrons with 
up and down spin states). This model has been studied 
by a number of different many-body techniques, such as 
self-consistent diagrammatic studies^, dynamical mean- 
field theory^ and renormalization group techniques^. The 
main message of all these studies was that this model 
contained the following basic physical properties: 

Upon decreasing the temperature below a certain char- 
acteristic T* , local electron pairing (not electron pair 
states!) begin to form and exist over a finite time as 
well as space interval. As a consequence of that, a pseu- 
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dogap centered around the Fermi level opens up in the 
density of states (DOS) of the itinerant electrons. Elec- 
trical transport is being controlled by relatively well de- 
fined single- electron states above T* , while below T* the 
electrons get heavily damped. Upon further decreasing the 
temperature, the electron pairs become better and better 
defined quasi-particles and ensure the transport, eventu- 
ally leading to a superconducting state controlled by phase 
fluctuations of the electron pairs. The resistivity, just be- 
fore entering this superconducting state, resembles that of 
an insulator, since it is determined by the electrons with 
their pseudogap features. 

In the present paper, we examine all these features 
within an exact diagonalization scheme, evaluating the 
static and dynamical properties of the BFM, by means of 
the Lanczos procedure extended to finite temperature^. 
This technique has been tested and successfully applied 
on different many-body systems of strongly correlated 
electrons^, and allows to determine the temperature de- 
pendence of various physical quantities, which, combined 
with dynamical properties, gives a complete information 
about the energy spectrum and excitations of the model 
under examination. In particular, within this scheme 
of computation the trace of thermodynamic expectation 
values is obtained by means of Monte Carlo sampling, 
while for the determination of the relevant matrix ele- 
ments the Lanczos method is used. 

We restrict ourselves to the symmetric case, where the 
bosonic level lies in the middle of the electronic band, i.e., 
As = 0. For a different position of the bosonic level the 
physics is qualitatively the same, the relevant parameter 
in this scenario being the concentration of the bosons. 
Due to their hard-core nature, they induce correlation 
effects which lead to qualitative changes in the physical 
properties as we go from low (ng > 0) to high (jib — 0.5) 
boson concentrations. 

Our local phase space being made up of eight states, 
a large computational memory is required. This limits 
the present study to maximally eight sites systems if the 
evaluation of the static and dynamical quantities does 
not involve transitions between subspaces with a differ- 
ent number of fermions and bosons. Moreover, for the 
study of the spectral functions for the fermions, together 
with the previous constraint we have to consider a situa- 
tion such that the discrete energy spectrum of our finite 
size system has contributions as close as possible to the 
Fermi level, in order to track the low energy features of 
interest here. Having the Fermi level lie in the center of 
the fermionic band (in order to avoid effects coming from 
any particular band structure), the maximal size of the 
system which we could consider is restricted to four sites 
rings allowing any concentration of fermions and bosons. 
As far as concerns the boson spectral function, we can 
investigate up to six sites rings since the constraint on 
the physical description of the fermionic properties does 
not apply. 

The relevant correlation functions determining the on- 
set of local electron pairing and spatial phase correla- 
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FIG. 1: Variation with temperature of the fermionic spectral 
function at the Fermi vector, evaluated for a 4-site ring with 
A B = 0, g/t = 0.2 and n tot = 8 (n B = 2). 



tions are given by (pf c il c il + h.c.) and (pf p~) and will 
be discussed in Section II. In Section III we examine the 
spectral properties of the hard-core bosons and compare 
them with the excitations of an equivalent XY model in 
the presence of a transverse field. The change-over from 
single-electron to electron pair transport will be discussed 
in Section IV on the basis of the temperature dependence 
of the optical conductivity and in section V in terms of 
persistent currents in mesoscopic rings. 



II. FROM LOCAL TO LONG-RANGE PHASE 
CORRELATIONS 



The spectral features of the fermions are determined 



by 
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where A F (k,uj) represent the electron and hole spectral 
functions, respectively, Q denotes the partition function 
and E m are the eigenvalues of the total Hamiltonian. 
In Fig. ^ we plot the evolution of the fermionic spectral 
function for electrons at the Fermi surface, i.e., k = kp. 



3 




FIG. 2: Temperature dependence of the on-site correlation 
function {pfrf), with t~ = Cij,Cjf, evaluated for a single site 
in the atomic limit (Fig. la). Comparison of this local corre- 
lation function (Fig. lb) and (p^ =0 p~ =0 ) (Fig. lc), measuring 
the long-range phase coherence of the hard-core bosons, eval- 
uated for an 8-site system for different boson concentrations 
riB. The value of the boson-fermion coupling is g/t — 0.4. 

We notice that with decreasing temperatures spectral 
weight is transferred from the Fermi energy to side bands 
which characterize the contributions from bonding and 
antibonding two-electron states. This emptying out of 
the spectral density near the Fermi energy is responsi- 
ble for the opening of a pseudogap in the DOS of the 
fermions below a certain characteristic temperature T* , 
which can be determined from the drop in spectral weight 
of Ap(kp, fi). 

The onset of local electron pairing, seen in the opening 
of a pseudogap in the electronic DOS, is manifest more- 
over in a strong increase of the local correlation func- 
tion between the hard-core bosons and the electron pairs 
(pf cyc,~f). This effect is already inherent in a single- 
site system and is illustrated in Fig. [SJi, showing a char- 
acteristic increase of this correlation function for tem- 
peratures of the order of <?, which is the energy gain 
of two electrons when they engage in a bonding state 
(1/ V2) [pf + 4 c + ] 1 0) . Going beyond the single-site case, 
this energy gain is reduced to g 2 /t as a consequence of 
the electron itinerancy entering into competition with the 
charge exchange mechanism. The intensity of the on- 
site correlation between bosons and electron pairs then 



correspondingly decreases, as illustrated in Fig. 2b for 
different total numbers of charge carriers, corresponding 
to ns = 1,2,3,4. Finally, in Fig. 2c we examine the 
long-range phase coherence for the hard-core bosons, de- 

scribed by <p+ =0 p- =0 ) - (l/N)YLMpiys)' We n0 " 
tice that this correlation function strengthens as T drops 
below a certain characteristic temperature (< T*). 

The evolution of T* and noticeably changes with 
the concentration of bosons. We derive these temper- 
atures from the inflexion point of these two correla- 
tion functions and present them as a function of ns in 
Fig. |3 While for small boson concentration T$ follows 
the same behavior as T*, it begins to show a down- 
wards trend upon further increasing the boson concen- 
tration. This confirms a recent result for the phase dia- 
gram of the BFM obtained within a diagrammatic tech- 
nique which explicitly takes into account the hard-core 
nature of the bosons^. In that study it was concluded 
that this non-monotonic behavior of should be related 
to the correlation effects in the hard-core boson subsys- 
tem, leading to an increasingly heavy mass of the bosons 
as their concentration approaches the dense limit, i.e., 
Ub = 0.5. This is to be expected for systems where 
the onset temperature of phase ordering is related to 
the phase stiffness ~ n p /m p of the electron pairs 
which is induced by the boson-fermion exchange cou- 
pling, n p = (nF^nFi) — (nF^)(nFi) denoting their density 
and m v their mass. 



III. BOSON SPECTRAL PROPERTIES: 
COMPARISON WITH XY PHASE 
FLUCTUATION SCENARIOS 



After having considered the spectral properties of the 
fcrmions in the preceding section, let us now discuss the 
spectral properties of the hard-core bosons. The appro- 
priate spectral function in this case is 



A B {q,uj) = A^(q,uj) - A B (q,u) 
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with = -7= J2i eZqTi pf ■ We ultimately want to com- 
pare this spectral function, evaluated within the BFM, 
with that of an effective XY model in a transverse field, 
described by 
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XY 



-'Ed 



Pi+6 



(4) 



The corresponding results are presented in Fig. show- 
ing the absolute values of these correlation functions (we 
recall that they are negative for negative frequencies tu). 
For the case considered here, a 6-site ring with 3 bosons 
and Ab = 0, we have a chemical potential /1 = and 
thus we choose h — in order to describe an equivalent 
physical situation within a XY model. The typical energy 
scale of the bosonic excitations is of the order of g 2 /2D 
which is the energy scale related to the exchange inter- 
action allowing the transfer of a boson from a given site 
to an adjacent one. The strong similarity between these 
two spectral functions over the entire range of q vectors 
is a strong indication that for low temperatures (T = 
in the present case) the excitation spectrum of the BFM 
is determined by phase fluctuations of the bosonic charge 
carriers. 

In Fig. we illustrate for the fully symmetric case 
(Ab = and n to t — 12 for a 6-site ring) how the spectral 
function of the hard-core bosons evolves as a function of 
temperature. We notice that upon increasing the tem- 
perature above T$, this spectral function becomes more 
and more incoherent. The quasi-particle features of the 
fcrmions and of the hard-core bosons develop in opposite 
directions as the temperature is lowered. Above T* the 
fermions are still rather well defined quasi-particles, while 
the hard-core bosons are completely incoherent. Between 
T* and the fermions lose their quasi-particle features, 
while those of the bosons are getting better and better 
defined. Finally, below the only good quasi-particles 
are the hard-core bosons and consequently also the local 
electron-pairs. This change in the quasi-particle features 
from single to two-particle states will become more ap- 
parent in the transport properties, discussed in the fol- 
lowing Section. 
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FIG. 4: Comparison at T = of the boson spectral func- 
tion Ab{q,oj) for ntot = 12 (ns = 3) on a 6-site ring with 
the spectral function of an equivalent XY model (denoted by 
AxY{q,u))- 



IV. FERMION VERSUS BOSON TRANSPORT 

The spectral properties of the fcrmionic single-particle 
excitations, examined in Section II, indicate a loss of 
spectral weight for wave vectors centered around kF, 
which is responsible for the opening of the pseudogap 
in the DOS. Full many-body calculations, involving self- 
energy corrections^, have shown that besides this loss 
of spectral weight, fermions also lose their quasi-particle 
features below T* . They become purely diffusive and 
show up in a resistivity which monotonically increases 
as the temperature tends to zero. At the same time, 
the bosons and respectively the electron pairs, which are 
diffusive above T*, become well defined quasi-particles 
upon approaching T^, as discussed in section III. Thus, 
the bosonic charge carriers with charge 2e are expected 
to ensure the transport upon approaching and even- 
tually lead to a drop in the resistance and ultimately to 
a superconducting state. 

These features can be illustrated within the present 
exact diagonalization study in terms of the optical con- 
ductivity arising from the fermions and the bosons, 
respectively. The electric field in the present model 
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FIG. 5: Evolution with temperature of the spectral properties 
of the hard-core bosons for a 6 site ring with 3 bosons (ntot = 
12) for two temperatures below and above T$. 

couples exclusively to the fermions, since the bosons 
are intrinsically localized. Nevertheless, because of the 
boson- fermion exchange coupling the bosons, and con- 
sequently also the fermion pairs, acquire itinerancy (as 
shown in Section III) and thus contribute to the trans- 
port via an Aslamazov-Larkin type contribution to the 
conductivity^ 3 . In order to capture this effect within the 
present exact diagonalization study, we decompose the 
current carried by the fermions jp — t £\ s a (cj ' a c i+ s a — 
h.c.) into a term associated with a current carried by the 

bosons jb = ts Yli s(Pi Pi+s ~ h-c.) plus a remainder jp 
such that 

j F {t)=3B{t)+lp{t). (5) 

Separating optimally the contributions coming from the 
current-current correlation function for the bosons from 
the remainder of the fermions, we require {j Bit) j Fit)) = 
0, which defines an effective boson hopping integral as 
ts = (jFit)jBit))/(j%it)). Such a decomposition is 
meaningful in a situation where the current transported 
by the fermionsjs rapidly dissipated as a function of time, 
described by (xfWxf^')); and its long time behavior 
is getting controlled by the current transported by the 
bosons, described by {j Bit) j Bit'))- Such a situation is 



realized for low temperatures, around and below T^. We 
can then safely neglect the contributions coming from the 
cross terms Of^j/b (£'))• 

We report in Fig. Et- the optical conductivity com- 
ing from exclusively the single-electron transport, given 
by (jFit)j F it')), and in Fig. |SJd the contribution com- 
ing from bosonic charge carriers. We notice that for the 
fermionic contribution ctf(w) to the optical conductivity 
we obtain a featureless metallic-type behavior for T > T* 
iT/t = 0.3 in our case), while upon going below T*, a 
pseudogap structure materializes. Upon further decreas- 
ing the temperature below T^, cf^) shows a well devel- 
oped gap, reflecting the onset of an insulating behavior, 
as it can be inferred from the temperature behavior of 
the dc conductivity, illustrated in the inset of Fig. [f)]i. 

At the same time, when sweeping through this same 
temperature interval, we notice the build-up of a struc- 
tured contribution in the bosonic component of the op- 
tical conductivity ctb(w), which develops precisely in the 
frequency interval where the pseudogap appears in the 
fermionic contribution o\p(w)- These features confirm the 
cross-over from fermionic to bosonic transport with de- 
creasing temperature. 

Of course, the spectral weight of the bosonic absorp- 
tion is much smaller with respect to that of the fermionic 
sector. The effective boson hopping is ts ~ g 2 /2D, 
and hence the relative spectral weight of the bosonic ab- 
sorption goes like ~ t 2 B . Nevertheless, it becomes clear 
from Fig. that below the transport at low frequency 
(ui/t < g 2 /2D) is completely due to bosons and pairs of 
fermions. 



V. PERSISTENT CURRENTS IN MESOSCOPIC 
RINGS 

The change-over between single-particle and two- 
particle transport is a salient feature of the BFM and 
its experimental verification would be vital in verifying 
such a scenario. As we have seen above, the onset of 
bosonic transport is not necessarily linked to a super- 
conducting state. It is sufficient to have phase-coherent 
free particle-like bosonic excitations. In mesoscopic rings 
where persistent currents are a characteristic feature^*, it 
should be possible to detect this change-over in transport 
when threading a magnetic flux through such a ring. This 
can give rise to persisting currents, periodically varying 
as a function of the flux $. In units of the flux quantum 
$0 = foc/e this periodicity should then be seen in multi- 
ples of $/$o = n if the transport is via single electrons, 
and of $/$o = n/2 if it is due to two-particle bosonic 
transport, n denoting an integer number. 

We investigate this problem in the standard way by 
evaluating the free energy F($>) in the presence of a 
magnetic field and then determine the resulting cur- 
rent as J/N = dF($)/d$. The behavior of fol- 
lows from the evaluation of the thermal average of the 
Hamiltonian (1), upon making the Peierls substitution 
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FIG. 6: Optical conductivity due to fermions (Fig. 4a) and 
due to bosons (Fig. 4b) for an 8-site ring with ntot = 16 
(jib = 4) and three characteristic temperatures, correspond- 
ing to above T* , between T* and and below T^. The 
displayed difference in the amplitude of the fermionic and 
bosonic conductivity is set by the ratio between the effective 
boson and fermion hopping (~ g 2 /2D). The scale is fixed by 
po — ha/e 2 where a is the atomic distance. 



t — > t cxp ijflf^J ■ We illustrate in Fig.0the variation 

in the ground state of F($) as a function of $/$o f° r 
the uncoupled system (g — 0) and the fully interacting 
system (g ^ 0). In the inset of the figure we also report 
the behavior of F(&) for small positive values of <f> at 
T = and at a finite temperature T > T$, in the case 
g/t = 0.4. As to be expected, for the case g = we find a 
situation where the applied field lowers the energy of the 
system, arriving at minima at $/$o = n /2 (see dotted 
line in Fig. [7J. The periodicity of one flux quantum of 
F($) suggests that the effect of the flux is to introduce a 
paramagnetic current made out of single-electron states. 
Upon switching on the boson-fermion interaction, g =/= 0, 
one finds that around $ = the free energy below 
increases with increasing $ before bending over and fol- 
lowing essentially the same behavior as that obtained for 
g = 0. This suggests that if a small flux is applied to the 
system, a diamagnetic screening current is building up 
which, eventually, gives way to a paramagnetic current 
stabilizing the system when the flux is further increased. 
Hence, upon switching on the boson-fermion interaction, 
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FIG. 7: Comparison of the free energy F(4>) as a function of 
the flux $ for the uncoupled (g = 0) and the coupled (g =fc 0) 
boson-fermion system. In the inset we present the evolution 
of F($) for positive $ close to 0, for T = and T = 0.04 
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FIG. 8: Evolution with temperature of the current induced 
by a magnetic field threading a mesoscopic ring, illustrating 
the change in periodicity of the current as we pass from tem- 
peratures above TU to below it. 
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the system evolves towards a state where the free en- 
ergy exhibits minima at values of the flux with periodic- 
ity $o/2- This is indicative of an induced current made 
out of carriers having charge 2e, and of an off-diagonal 
long-range order in the ground state. These conclusions 
are based on results derived a long time agoi^^ in or- 
der to provide a criterion to distinguish between normal 
and superconducting ground states from the functional 
form of i r ($), without any knowledge of the symmetry of 
pair-pair correlation functions. The fact that the mini- 
mum of the free energy at <&/$() — is not the same as at 
$/$o = 1/2 is due to the finite-size calculation presented 
here, analogous to the study of the anomalous flux quan- 
tization in the negative U Hubbard problem (see Fig. 1 
in reftf). As the temperature is increased above T^, we 
notice that the minimum of F at $ = gradually disap- 
pears, thus stabilizing the paramagnetic single-electron 
current. 

In Fig. [SJwe present the current induced by the mag- 
netic flux threading an 8-site mesoscopic ring as a func- 
tion of ^/^oj f° r three characteristic temperatures, two 
below and one above T^, We again can see the change- 
over from two-particle to single-particle transport. Ex- 
periments on mesoscopic ring made out of a high-T c com- 
pound might possibly have seen such two-particle trans- 
port at low temperatures^. 

VI. CONCLUSIONS 

The present exact diagonalization study of the BFM 
confirms previous results obtained by various many-body 
techniques, all of which have their built-in shortcomings. 
The competitive effect of local and long-range correla- 
tions is a new element in the present approach, allowing 



to extract some relevant physics underlying the BFM . 
They are characterized by the energy scales ksT* and 
ksTfi, describing respectively the onset of local pairing 
and the establishment of a phase-correlated state on a 
long range, the latter being indicative of a superconduct- 
ing phase in an infinite system. The change-over from 
single-particle to two-particle transport is most clearly 
seen and discussed in terms of the optical conductivity, 
showing below T* a single-particle pseudogap feature to- 
gether with a two-particle contribution developing in the 
frequency window opened up by the pseudogap. A clear- 
cut experiment for verifying the present scenario would 
be to study the change-over form single-electron to two- 
electron transport in mesoscopic rings, the theoretical 
expectations for which have been discussed. 

In the present study, we limited the analysis to the case 
of optimal resonant scattering between the electrons and 
the local bound electron pairs, this being the physically 
most interesting case for the interplay between phase and 
amplitude fluctuations. Such a situation is realized when 
the bosonic level Ab lies in the middle of the fermionic 
band (i.e. As = 0). By changing the position of the 
bosonic level, there are no qualitative changes as far as 
concerns the static and dynamical response and the ap- 
pearance of two characteristic energy scales which control 
the short and long range phase coherence. In particular, 
the three- peak structure in the density of states (Fig. 1) 
remains the same, but the spectral weight is now asym- 
metric with respect to the chemical potential. 
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